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FOREWORD 


This  interim  report  presents  details  of  the  first  tunnel  entry  of  an 
atmospheric  probe  model.  The  project  is  a  joint  effort  of  the  Air  Force 
Cambridge  Research  Laboratories  (AFCRL),  Arnold  Engineering 
Development  Center  (AEDC),  and  Mithras,  Inc. ,  Cambridge,  Massa¬ 
chusetts.  The  investigation  is  being  conducted  by  D.  Golomb  and 
F.  P.  DelGreco  of  AFCRL,  J.  A.  F.  Hill  and  R.  E.  Good  of  Mithras, 
Inc.,  and  J.  A.  van  der  Bliek  and  R.  A.  Cassanova  of  AEDC  under 
Program  Element  62405424,  Project  7635. 

The  results  of  tests  presented  were  obtained  by  ARO,  Inc.  (a  sub¬ 
sidiary  of  Sverdrup  &  Parcel  and  Associates,  Inc.  ),  contract  operator 
of  the  AEDC,  Air  Force  Systems  Command  (AFSC),  Arnold  Air  Force 
Station,  Tennessee,  under  Contract  AF  40(600)- 1200.  The  tests  were 
conducted  from  September  1  to  23,  1965,  under  ARO  Project  No.  SB0512, 
and  the  manuscript  was  submitted  for  publication  on  April  25,  1966. 

Final  results  will  be  published  under  joint  authorship  of  the  above 
investigators. 

The  numerical  computations  of  the  radial  inversion  described  in 
Appendix  VI  were  carried  out  by  W.  C.  Armstrong  of  Central  Computing 
Operations,  ARO,  Inc. 

This  technical  report  has  been  reviewed  and  is  approved. 

James  N.  McCready  Leonard  T.  Glaser 

Major,  USAF  Colonel,  USAF 

AF  Representative,  AEF  Director  of  Test 

Directorate  of  Test 
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ABSTRACT 


Low  density  wind  tunnel  tests  were  carried  out  at  Mach  3  to 
simulate  the  release  of  nitric  oxide  in  the  upper  atmosphere  from 
a  sounding  rocket.  The  supersonic  flow  contained  atomic  oxygen, 
and  the  chemiluminescent  reaction  of  O  with  NO  was  investigated. 
This  interim  report  presents  some  of  the  details  of  the  experimental 
arrangement  and  preliminary  results  of  the  first  test  series. 
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SECTION  I 
INTRODUCTION 


In  connection  with  the  AFCRL  upper  atmosphere  research  program, 
attempts  have  been  made  to  determine  the  concentration  of  atomic 
oxygen  (O)  in  the  upper  atmosphere  (90  to  140  km  altitude).  One  method 
(Ref.  1)  consisted  of  releasing  nitric  oxide  (NO),  at  a  known  rate,  from 
a  rocket  at  night.  The  light  from  the  trail  was  analyzed  using  ground- 
based  cameras.  The  method  is  based  on  the  overall  reaction: 

NO  +  0  — - *  N02  +  hv  (1) 

If  the  reaction  rate  constant  is  known  from  laboratory  experiments, 
the  concentration  (0]  can  be  determined  if  the  light  output  and  are 

known.  It  is  suggested  in  Ref.  2  that  the  simple  direct  reaction  of 
Eq.  (1)  has  to  be  replaced  by  a  scheme  of  several  reactions  involving 
various  intermediate  excited  states.  In  fact  the  spectrum  associated 
with  the  reaction  is  almost  a  continuum  (Fig.  1),  suggesting  a  very 
large  number  of  excited  states.  Using  the  experimental  rate  constant 

K,  -  6.4  x  10""  - cm? - 

molecule  -  sec 

from  Ref.  3  and  the  known  I  and  -  ,  the  O  concentration  was  deter- 

dt 

mined  in  Ref.  1  from  three  rocket  flights.  Thus: 

1  =  =  K‘  £  tN0]  t°]  ( 2) 

where  £  is  a  flow  field  function.  It  was  then  found  that  typically  at 
120  km  altitude,  the  O  density  was  2500  times. the  total  atmospheric 
density.  Ref.  1. 

The  purpose  of  the  present  project  is  to  simulate  free-flight  condi¬ 
tions  with  the  aim  of  examining  the  flow  structure  around  the  rocket  and 
to  measure  the  light  output  in  a  stream  with  a  known  O  concentration  for 
a  given  rate  of  NO  release.  It  is  hoped  that  this  will  allow  a  more 
realistic,  quantitative  interpretation  of  the  rocket  flight  data. 

This  report  concerns  the  experimental  arrangement  and  some 
preliminary  results  obtained  in  August  and  September,  1965.  A 
second  test  series  is  planned  in  the  spring  of  1966  to  complete  this 
project. 
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SECTION  II 

SIMULATION  OF  UPPER  ATMOSPHERE  ROCKET  FLIGHT 


In  order  to  produce  meaningful  test  results ,  the  flow  pattern  around 
the  model  should  be  aerodynamically  and  chemically  similar.  The 
radius  of  the  headglow  (Section  VII)  produced  by  the  NO  injected  in  the 
upstream  direction  was  taken  as  the  typical  dimension.  For  the  vehicle 
"Dinah",  Ref.  1,  this  radius  (Rc)  is  tabulated  below  as  a  function  of 
altitude: 


Z,  km 

100 

120 

140 

Rc,  m 

11 

22 

40 

Pw,  kg/m8 

4. 97  x  10'7 

2.  44  x  10"8 

3. 3  x  10"9 

A„,,  m 

0.  16 

3.  23 

22.  5 

471 

24.  7 

2.  89 

M* 

4.  1 

2.  5 

1.  3 

K 

220 

440 

800 

Some  of  the  other  free-flight  parameters  are  also  listed. 

In  selecting  the  tunnel  conditions,  it  is  desirable  to  duplicate  the 
Mach  number  and  Reynolds  number.  This  includes  duplication  of  the 
Knuds en  number 


K  =  1.485  J)L_  for  y  =  1.4 

n  Re 


The  radius  of  the  headglow  in  the  test  section  is  typically  (Rc)  ~  5 

(R  ),, 

and  this  fixes  the  scale  factor  K  =  -  —  listed  above.  Since,  for  aero- 

(Rc>wt 


dynamic  simulation  (Re){{  =  (Re)wt  and  for  chemical  simulation  it  is 
desirable  to  produce  the  same  temperature,  the  density  is  scaled  by  a 
factor  K.  Thus  the  product  Rc  is  identical  for  wind  tunnel  and  free- 
flight  conditions.  If  the  dominant  chemical  reactions  are  bimolecular, 
the  flow  will  be  chemically  similar  as  was  shown  for  air  in  Ref.  4.  * 


*It  is  interesting  to  note  that  binary  scaling  cannot  be  generally 
applied  in  wind  tunnel  testing  except  perhaps  in  some  free-flight  ranges. 
The  reason  is  that  models  have  to  be  scaled  down  and  the  density  in  the 
test  section  cannot  be  increased  accordingly. 
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For  the  above  tabulated  data,  the  desired  tunnel  conditions  are 
plotted  as  a  function  of  simulated  altitude  in  Fig.  2.  The  regime  of 
duplication  of  each  of  these  parameters  in  the  present  chamber  using 
the  Mach  3  nozzle  and  with 

T0  =  290  to  800° K 
Po  =  50  to  400  y.  Hg 

is  shown  by  the  double  lines  in  Fig.  2. 

In  the  wind  tunnel  tests  a  mixture  of  N2  and  O2  flowed  through  the 
nozzle.  The  oxygen  was  partially  dissociated. 

The  gas  flow  from  the  sounding  rocket  consisted  of  pure  NO  but 
during  the  wind  tunnel  tests  a  mixture  of  NO  and  N2  was  used.  By  keep¬ 
ing  the  total  model  pressure  constant  and  varying  the  NO  content,  the 
headglow  geometry  remained  constant  but  the  concentration  of  NO  in  the 
reaction  regime  was  varied  from  100  percent  down  to  a  fraction  of  a 
percent. 


SECTION  III 

PRODUCTION  OF  ATOMIC  OXYGEN 


Several  methods  of  producing  O  in  the  supersonic  stream  were 
considered. 

3.1  PHOTO  DISSOCIATION 

Dissociation  of  O2  takes  place  in  the  upper  atmosphere  by  absorp¬ 
tion  in  the  ultraviolet.  There  is  a  broad  continuum  (Schumann-Runge 
band)  with  a  maximum  near  0.  1450  m.  At  lower  wavelengths,  there 
are  several  strong  absorption  peaks  (Ref.  5).  Since  this  method  would 
tend  to  produce  O  in  the  same  state  as  in  the  upper  atmosphere,  it 
would  seem  worthwhile  to  explore  the  possibility  of  using  ultraviolet 
radiation  in  a  wind  tunnel. 

3.2  ARC  HEATER 

Arc  heaters  would  produce  a  gas  with  dissociated  O2.  It  was 
thought,  however,  that  without  considerable  development  a  low  density 
arc  jet  would  contain  considerable  contamination  from  the  electrodes 
and,  since  the  arc  itself  would  operate  at  a  much  higher  temperature 
than  the  desired  temperature  (~1000°K),  other  gas  species  would  be 
present.  It  is  noted  that  the  percentage  of  O  at  equilibrium  conditions 
of  air  at  1000°K  is  negligibly  small  so  that  the  O  content  would  have  to 
be  frozen  from  a  higher  temperature  condition. 
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3.3  RF  GENERATOR 

Radio  frequency  excitation  has  been  used  to  dissociate  O2  (Ref.  6). 
Laboratory  experiments  at  AFCRL  and  also  at  AEDC  indicated  con¬ 
siderable  interference  with  nearby  instrumentation.  In  a  separate 
laboratory  arrangement,  similar  to  that  used  here,  an  rf  generator 
operating  at  10  me  was  tried.  It  produced  nominally  a  250  vrms  signal 
in  a  wire  wrapped  (7  to  13  turns)  around  a  13-mm  glass  tube. 

The  rf  generator  set  produced  a  color  more  toward  red  (almost 
pink)  than  the  radar  set  (see  below).  The  luminous  region  extended  far 
downstream  of  the  discharge  region.  The  reaction  tube  could  be  filled 
almost  completely  with  a  soft  red  to  pink  afterglow.  No  meaningful 
measurements  could  be  made  since  the  electromagnetic  field  affected 
all  electrical  instrumentation.  No  attempt  was  made  to  shield  the  in¬ 
strumentation. 


3.4  MICROWAVE  GENERATOR 

The  method  of  dissociation  of  O2  with  microwaves  has  been  used  in 
several  chemical  laboratories  (Ref.  7).  Pilot  experiments  were  carried 
out  for  this  test  at  AFCRL  and  AEDC.  The  microwave  unit  used  at 
AEDC  consisted  of  a  converted  pulsed  radar  set.  The  waveguide  termi¬ 
nated  in  a  horn  antenna  shown  in  Fig.  3.  The  set  has  the  following  char¬ 
acteristics: 

Average  Power  =  60  w 
Triggering  Frequency  =  1000  cps 
Peak  Power  ®  300,  000  w 
Output  Frequency  =  9.7  kmc  (X-band) 

With  only  O2  flowing  through  a  9-mm  glass  tube,  the  degree  of  dissocia¬ 
tion  was  from  0.  2  to  0.  3  percent  at  mass  flow  rates  from  3  to  200  mg /sec. 
In  the  AFCRL  experiments  a  continuous -wave  microwave  generator  oper¬ 
ating  at  2.  5  kmc  and  feeding  into  a  resonant  cavity  around  a  quartz  tube 
was  used.  In  view  of  the  success  and  ease  of  operation  of  this  type  of 
microwave  generator,  this  arrangement  was  used  in  the  wind  tunnel 
experiments. 


SECTION  IV 

LOW  DENSITY  WIND  TUNNEL 


The  tests  were  conducted  in  the  Aerospace  Research  Chamber 
(ARC)  (8V)  of  the  Aerospace  Environmental  Facility  at  AEDC.  The 
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chamber,  operated  as  a  low  density  wind  tunnel,  is  shown  in  Fig.  4. 

The  chamber  is  8  ft  in  diameter  and  16  ft  long.  The  nozzle  is  mounted 
on  a  removable  dished  head.  The  pumping  equipment  of  the  chamber 
consists  of  a  6 -in.  oil  diffusion  pump  in  series  with  a  60-liter /sec 
mechanical  pump. 

The  major  pumping  action  during  a  tunnel  run  is  provided  by  a 
series  of  radially  positioned  cryopanels  cooled  with  gaseous  helium 
(GHe)  at  from  15  to  20°K  with  a  total  area  of  240  ft^.  The  GHe  is  cooled 
in  a  continuously  operating  1-kw  refrigerator.  The  cylindrical,  per¬ 
forated  precooler  and  the  chamber  end  panel  are  cooled  by  liquid 
nitrogen  (LN2)  at  77°K.  The  cylindrical  part  of  the  chamber  is  pro¬ 
vided  with  an  LN2  liner  which  serves  as  a  heat  shield  for  the  helium 
cryopump. 

The  nozzle,  shown  in  Fig.  4,  is  cooled  with  LN2  to  reduce  the  noz¬ 
zle  wall  boundary-layer  thickness.  During  this  test  a  Mach  3  nozzle 
was  used. 

During  the  present  test,  O2  and  N2  were  fed  to  the  gas  heater  from 
standard  bottles,  via  standard  flowmeters.  The  gas  is  passed  through 
two  separate  tubular  gas  heaters  into  a  heated  plenum  chamber.  The 
O2  was  passed  through  a  quartz  manifold  with  five  microwave  cavities 
to  dissociate  part  of  it.  When  these  cavities  were  powered  by  five 
80-w  microwave  generators  (2.  5  kmc)  the  O2  heater  upstream  of  the 
cavities  was  not  turned  on.  During  operation,  the  cavities  became  hot 
because  of  electric  power  dissipation  and  had  to  be  cooled  with  com¬ 
pressed  air.  A  closeup  of  this  manifold  is  shown  in  Fig.  3b.  The 
cavity  size  and  thus  the  quartz  tube  size  is  fixed  by  the  frequency  em¬ 
ployed.  Furthermore,  the  operating  pressure  in  the  cavities  appeared 
to  be  limited  to  10  mm  Hg.  This  led  to  the  requirement  of  five  parallel 
cavities  with  five  power  supplies  instead  of  one  larger  power  unit.  The 
quartz  tube  downstream  of  the  cavities  was  lined  with  Teflon®  to  reduce 
O  recombination. 

The  plenum  chamber  temperature  was  held  at  the  desired  reser¬ 
voir  temperature.  It  is  noted  that  the  pressure  in  the  tubular  gas  heater 
coils  was  from  1  to  10  atm  to  provide  effective  gas  heating  which  could 
not  be  accomplished  in  the  plenum  chamber  operating  at  from  50  to 
500  p  Hg. 

Several  control  valves  were  installed  to  achieve  the  desired  pres¬ 
sures  at  various  positions  in  the  gas  inbleed  system. 
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The  operating  range  during  the  present  test  was 

Pc  =  50  to  400  fi  Hg 


To  =  290  to  800°  K 


pM=  1  to  5  n  Hg 
=  85  to  27 5° K 
Re^/cra  =  1.2  to  33 
=  3.2  to  3.5 
A  =  0.2  to  3  cm 

oo 

Uniform  Core  Diameter  =  10  to  30  cm 


Extensive  calibration  results  will  be  published  in  a  future  report.  The 
test  conditions  are  listed  with  each  run  reported  in  this  paper. 


SECTION  V 

MEASUREMENT  OF  ATOMIC  OXYGEN  CONCENTRATION 


The  percentage  of  O  produced  by  the  microwaves  can  be  deter¬ 
mined  by  means  of  a  titration  technique,  developed  by  Kaufman  (Ref.  7). 

A  pilot  apparatus  was  built  with  the  object  of  gaining  experience  with 
this  technique,  which  was  almost  identical  to  that  shown  in  Fig.  5. 

Oxygen  was  admitted  through  a  regulator  valve  into  a  small  glass  tube 
where  dissociation  took  place.  From  there  the  gas  entered  a  reaction 
tube  in  which  either  NO2  or  NO  was  admitted  from  a  bottle.  The  reac¬ 
tion  of  interest  then  took  place  in  this  large  tube,  and  the  gases  were 
pumped  to  the  atmosphere  from  this  tube  by  means  of  a  mechanical  pump. 

The  velocity  and  pressure  in  the  reaction  tube  were  regulated  by 
means  of  a  regulator  valve  after  the  flowmeter  and  by  the  valve  before 
the  pump. 

Pressures  and  flow  rates  were  measured  at  points  shown  in  Fig.  5. 

The  mass  flow  rates  of  NO  and  NO2  were  determined  from  the  re¬ 
corded  pressure  decay  of  a  5 -liter  flask  which  was  filled  to  a  pressure 
up  to  20  mm  Hg  before  a  run.  The  NO2  was  obtained  from  a  bottle  of 
liquid  N2O4.  The  liquid  evaporates  and  dissociates  below  atmospheric 
pressure.  The  equilibrium  composition  of  the  system, 

2N02  t=;  N204 
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is  considered  in  Appendix  I.  In  deducing  the  mass  flow  rate  of  NO2, 
the  degree  of  dissociation  was  taken  into  account;  the  appropriate 
formulas  are  derived  in  Appendix  II. 

The  following  simplified  equations  describe  the  reaction  that  occurs 
when  NO2  is  admitted  to  the  reaction  tube  containing  O2  and  O: 

0  +  N02  — *  NO  +  02  (fast)  ^2) 

NO  +  0  —  N02  +  hv  (slow)  ^ 

It  was  shown  in  Ref.  7  that  the  presence  of  third  bodies  such  as  O2 
and  N2  did  not  affect  the  reaction  rate  and  the  light  output. 

As  the  NO2  flow  rate  is  increased,  fewer  O  atoms  are  available  for 

the  light  producing  reaction  of  Eq.  (1)  and  finally  when  d  ^  =  d 

the  reaction  of  Eq.  (1)  cannot  take  place.  At  this  titration  point,  the 
light  is  cut  off  sharply.  Visual  observation  of  the  cutoff  point  appeared 
to  give  the  same  results  as  the  photomultiplier  for  the  present  arrange¬ 
ment. 

It  is  of  interest  to  note  that  the  same  technique  can  be  used  to  deter¬ 
mine  the  concentration  of  N  in  a  flow  (Refs.  8  and  9).  In  this  case,  NO 


is  injected  into  the  flow: 

N  +  NO  — »  0  +  N2  (4) 

N  +  0  +  M  — >  NO  +  M  +  hv  (5) 

NO  +  0  — »  N02  +  hv  (1) 

NO  +  0  +  M  — ♦  N02  +  M  (6) 


Here  Eq.  (4)  is  much  faster  than  the  other  reactions  so  that  the  N  flow 
produces  an  equivalent  O  flow,  resulting  in  a  light  output  because  of  its 
reaction  with  NO.  At  the  titration  point,  the  N  and  NO  are  both  fully 
consumed  before  the  reactions  of  Eqs.  (5),  (1),  or  (6)  can  take  place, 
resulting  in  zero  light  output. 

Although  the  method  worked  well  in  the  pilot  apparatus  (Fig.  5),  it 
was  not  possible  to  extinguish  the  light  completely  when  this  reaction 
tube  was  mounted  in  the  wind  tunnel  (Fig.  6).  This  was  thought  to  be 
caused  by  the  higher  velocity  and  lower  pressure  in  the  tube  when 
mounted  in  the  wind  tunnel.  The  velocity  in  the  tube  is  determined 
partly  by  the  intake  area.  Several  ratios  of  intake  area  versus  tube 
area  were  tried.  When  this  ratio  was  smaller  than  0.  1  the  mass  flow 
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through  the  tube  became  so  small  that  the  light  output  became  very  weak. 
The  calculation  of  the  mass  flow  through  the  tube  and  the  intake  effi¬ 
ciency  of  the  final  configuration  are  discussed  in  Appendix  III. 

It  was  noted  that  when  NO  was  admitted  to  the  sampling  tube,  the 
light  output  was  sufficient  in  all  cases  to  be  recorded.  Therefore  the 
following  method  of  calibrating  the  tunnel  flow  for  local  O  concentration 
was  developed. 

The  sampling  tube  was  first  calibrated  outside  the  tunnel  in  the 
arrangement  of  Fig.  5.  For  a  given  flow  rate  of  O2  and  power  setting 
of  the  microwave  generator,  the  concentration  of  O  was  determined 
with  the  NO2  titration  technique.  Then,  for  the  same  settings,  NO 
instead  of  NO2  was  admitted  through  the  Teflon  loop  (Fig.  5).  The 
light  output  of  the  reaction  of  Eq.  (1)  was  recorded  with  a  photomulti¬ 
plier  tube  (PM).  In  this  manner,  the  light  output  was  obtained  for 
various  degrees  of  O2  dissociation.  This  constituted  the  calibration  of 
the  sampling  tube.  The  tube  was  next  mounted  in  the  wind  tunnel  and 
from  the  measured  PM  reading  the  local  degree  of  O2  dissociation  in 
the  test  section  was  determined.  When  this  tube  was  mounted  in  the 
tunnel,  N2  from  the  nozzle  flow  was  also  present  in  the  tube,  whereas 
the  sampling  tube  calibration  was  carried  out  with  O2  only.  To  account 
for  this  and  the  added  NO  mass  flow  in  the  sampling  tube,  a  correlation 
formula,  derived  below,  was  used. 


The  method  is  based  on  two  assumptions: 

1.  The  concentration  of  O  in  the  reaction  zone,  in  particular  at 
the  PM  station,  is  the  same  as  that  in  the  free  stream.  The  effects  of 
wall  and  volume  recombination  are  discussed  in  Appendix  IV  for  an  ideal, 
constant  velocity,  constant  pressure,  glass  tube,  and  it  is  shown  that  for 
the  present  configuration  (u  >  10  m/sec  and  P  <  0.  5  mm  Hg)  the  recom¬ 
bination  of  O2  is  expected  to  be  less  than  a  few  percent. 

2.  The  reaction 

NO  +  0  - *  N02  +  photon 

is  so  slow  that  only  a  negligibly  small  amount  of  O  is  consumed  by  this 
reaction  at  the  PM  station.  Considering  the  rate  of  [0]  consumption 


d  [0] 

dt 


=  K,  [NO]  [0] 


and  assuming  [NO]  »  [0]  and  therefore  [NO]  =  constant,  as  was  the  case 
in  the  present  tests,  the  concentration  of  O  is  given  by 

-[01-  „  exp/-K,  [NO] 

[0],  l  u  J 

in 
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The  requirement  is  therefore 

Ki  [NO]  x_  «  1  or  JL  »  K,  [NO] 

U  X 

Typically  ^-  <  10"2  so  that, 

assuming  K.  =  6.4  x  10  ~17 - - -  (Ref.  3),  this  gives 

particle-see 

[NO]  «  1.6  x  10ls  molecules/cm 3 
or 

P  =  .MNQ  [NO]  «  8  x  10" 5  gm/cm3 
N  0  ^ 

At  room  temperature  this  corresponds  to  a  pressure  of 

PN0  =  Pno  yp-  T  «  7  x  io4  dynes/cm2  -  5  x  io4  ^  Hg  which  was  always 

M  N  0 

satisfied. 

From  Eq.  (1)  the  photon  production  is 

I  =  =  Kj  [NO]  [0]  (?) 

at 


(8) 


(9) 


I 
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For  NO 


But  also. 

Assuming  uNQ  =  utot  ,  this  gives 


NO  “  (PNO  uNO  4  d2)st 
mst  =  hot  “tot  ~  d2 


NO 


“NO  _  p  “*st 


tot 


or 


rfi 


Pun  P 


NO 


(10) 


NO  '  tot  m 

s  t 

In  the  sampling  tube,  say  at  the  PM  station,  the  density  is 

Plot  =  P  o2  +  Pn2  +  Po  +  Pno 

In  all  cases,  Pq  is  very  small  compared  to  p^,  but  P]^o  ^  noi;  always 
negligible  so  that 


P o2  +  Pn2  -  Plot  -  Pno 


1  - 


NO 


tot 


Substitution  of  the  last  expression  in  Eq.  (9)  and  then  in  Eq.  (8),  together 
with  Eq.  (10)  gives 


I  =  Kx 


mnomo 


ai 3 


and  the  correlation  formula  is 

I  KN2 


(11) 


The  sampling  tube  temperature,  and  therefore  the  gas  temperature 
as  can  be  shown  easily,  was  kept  the  same  during  calibration  and  tunnel 
runs  (heating  tape  controlled  by  a  thermostat  was  wrapped  around  the 
tube  as  shown  in  Fig.  6).  If  p^ot  were  replaced  by  Ptot>  which  was 
measured,  the  maximum  error  going  from  pure  O2  to  pure  N2  would  be 


1.  31.  Therefore,  the  difference  in  molecular  weight  has  to  be 


taken  into  account. 


This  is  done  by  writing 


PM 

P  Ptot 


P 


to  t  M,tf 

RT 


10 


AEDC-TR-66-105 


The  final  correlation  formula  is  now,  for  a  given  temperature: 


NO 


1  - 


NO  T}2 


PM  Ueff 


(12) 


The  pressure  at  the  PM  location,  PpM>  was  determined  from 


P2 


P  M 


P,2 


256 


m 


St 


A  x 

Md4 


as  follows  directly  from  Eq.  (Ill - 1 )  in  Appendix  III.  Finally  it  is  noted 
that  the  correct  procedure  to  calculate  would  be  to  start  from 
Dalton's  law: 


Pi 


P2 


Pi 

m2 


RT 


P. 


RT 


where  Piot  =  1  P.  .  This  would  give  for  the  density 

£  p  m 
p.  .  _  i _ i  i 


In  the  sampling  tube,  the  values  of  Pi  are  not  known.  In  the  present 
tests,  it  was  assumed  that 

f  P,  M,  =  p.ot  M.ll 

where 

mo  +  iil,  Mn 

Meff  =  -  2  ■- - — -  =  ^3  M0  +  (1-/3)  Mn  =  4(7  +  (3) 

m02  +  mN2  2 


The  results  of  the  sampling  tube  calibration  are  shown  in  Fig.  7. 
This  figure  shows  that  the  PM  output  was  proportional  to  the  degree  of 
dissociation  as  indicated  by  Eq.  (12). 

The  results  of  a  series  of  calibrations  in  the  tunnel  in  the  configura¬ 
tion  shown  in  Fig.  6  are  given  in  Fig.  8.  For  each  run  the  value  of  T, 
from  Eq.  (12),  was  calculated  and  then,  using  Fig.  7,  the  value  of  a  was 
determined. 

The  value  of  a  decreases  with  increasing  mQ2  for  a  given  reservoir 
pressure  (Fig.  8b).  Apparently  the  net  input  of  dissociation  energy  is 
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about  the  same  and  so,  when  the  mass  flow  through  the  microwave 
cavities  is  increased,  the  degree  of  dissociation  decreases.  For  a 
given  mass  flow  through  the  cavities,  the  degree  of  dissociation  in  the 
test  section  increases  with  increasing  PQ  (Fig.  8a).  A  possible  ex¬ 
planation  may  be  found  in  the  geometry  of  the  inbleed  system. 


The  quartz  inbleed  line  is  blocked  off  at  the  end  and  the  gas  (02 
and  O)  enters  through  a  series  of  holes  on  the  circumference  of  the 
cylindrical  tube.  Presumably  the  O  spreads  more  at  low  N2  flow  rate, 
and  more  recombination  takes  place  at  the  heated  settling  chamber  walls. 
During  calibration  runs,  the  sampling  tube  was  moved  laterally  several 
inches  which  produced  no  change  in  the  PM  reading.  It  is  therefore 
assumed  that  the  O  concentration  was  uniform  in  the  test  area. 

For  a  given  plenum  pressure  PQ  the  degree  of  dissociation  is  much 
lower  at  T0  =  800°K  than  at  T0  =  290°K.  This  is  partly  because  a  much 
lower  Nj>  flow  rate  is  required  for  a  given  m.02  to  achieve  the  same  PQ 
at  a  higher  temperature. 

Finally  the  microwave  cavities  were  turned  off  one  by  one,  and  the 
degree  of  dissociation  was  proportional  to  the  number  of  cavities  in 
operation  (Fig.  8c),  showing  that  each  cavity  produced  about  the  same 
amount  of  O. 


If  D  is  the  energy  of  dissociation  per  molecule,  then  the  total  dis¬ 
sociation  power,  Pdiss>  required  is 


Pdiss  =  «  M—  ND-  W 


Using  D  =  5.  080  ev  =  8.  14  x  10“19  joules /molecule,  and  assuming  uni¬ 
form  dissociation  at  the  nozzle  exit  plane,  Pdiss  =  40  to  80  w  for  the 
present  run  series.  The  cavities  were  set  at  90  percent  of  full  power. 
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which  gives  360  w  output  from  the  microwave  generators.  The  overall 
efficiency  is  therefore  from  10  to  20  percent.  Using  100-percent  power 
(400  w)  the  maximum  degree  of  dissociation  for  mQ2  =  0.  21  gm/sec 
would  be  12.  4  percent. 


SECTION  VI 

MODELS  AND  EXPERIMENTAL  ARRANGEMENT 


Two  models,  A  and  B,  were  used  to  simulate  the  rocket  nose  from 
which  NO  is  released  into  the  atmosphere.  The  dimensions  of  the  models 
are  shown  in  Fig.  9.  Model  A  which  was  similar  to  the  rocket  test 
vehicles  was  provided  with  orifices  on  the  second  conical  surface,  and 
model  B  had  an  orifice  at  the  apex  of  the  cone.  Rockets  of  the  model  B 
type  have  not  yet  been  flown.  The  models  were  connected  via  a  supply 
line  to  gas  bottles. with  a  mixture  of  NO  and  N2.  The  pressure  inside 
the  model,  measured  with  a  strain-gage  pressure  transducer,  was  regu¬ 
lated  with  a  regulator  valve  near  the  bottles.  Model  pressure  (PDj) 

ranged  from  1  to  20  atm.  By  varying  the  mixture  of  NO  and  N2  in  the 
model  the  same  headglow  geometry  (P0j  and  POc0  constant)  could  be  main¬ 
tained  while  varying  the  concentration  of  NO  in  the  reaction  region. 

The  overall  light  output  of  the  headglow  was  measured  with  a  PM 
mounted  on  the  tunnel  sidewall. 

The  headglow  was  photographed  with  35 -mm  Kodak  Tri-X  film  with 
an  exposure  time  of  0.  1  sec  and  an  f/2  lens  opening.  Color  motion 
pictures  were  taken  during  some  of  the  runs  using  high-speed  Ekta- 
chrome®  film  with  1/6 -sec  exp  and  f/0.  95. 

The  film  and  PM  were  calibrated  with  a  light  box,  which  was  mounted 
at  the  position  of  the  headglow.  This  light  box,  in  turn,  was  calibrated 
against  a  standard  source  at  AFCRL. 


SECTION  VII 

EXPERIMENTAL  RESULTS 

7.1  CORRELATION  LIGHT  OUTPUT -MODEL  A 

Typical  photographs  of  the  headglow  with  Model  A  are  shown  in 
Fig.  10.  It  is  clear  that  the  mixing  and  reaction  zone  is  confined  to  a 
curved,  thin,  sheet -like  volume.  The  light  output  was  measured  from 
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the  side,  but  presumably  the  PM  reading  would  depend  on  its  azimuthal 
position  with  respect  to  the  model. 

The  results  of  the  runs  with  Model  A  are  tabulated  in  Table  I.  All 
runs  were  made  with  a  tunnel  reservoir  temperature  of  800°K  and  a 
model  (NO)  reservoir  temperature  of  290°K.  Only  pure  NO  was  re¬ 
leased  from  the  model.  The  purpose  was  to  simulate  only  previous 
rocket  flights. 

The  light  output  is  plotted  in  Fig.  11a  for  constant  tunnel  reservoir 
conditions  as  a  function  of  the  NO  reservoir  pressure,  P0j.  Several 

attempts  were  made  to  plot  these  data  in  a  unified  manner  (Figs,  lib,  c, 
and  d).  As  a  first  approximation,  the  geometry  of  the  headglow  is 
primarily  a  function  of  P04/P0  .  Plotting  the  light  output  versus 

Po-j/Po*,  gave  no  systematic  difference  for  curves  of  different  (Pcducd)q 

mass  flow  of  O  per  unit  area,  as  shown  in  Fig.  lib.  The  variation  of 
( p ^ u^o )q  was  small  and  the  points  were  all  within  25  percent  from  the 
mean  line.  Indeed,  plotting  of  light  output  divided  by  O  mass  flux 
against  Poj/Po®  gave  essentially  the  same  spread  of  data  points 

(Fig.  11c).  In  Fig.  lid,  the  same  parameter  was  plotted  as  a  function 
of  Po-i/P®*  but  this  gave  a  slightly  larger  spread  of  points. 

J 

Finally  the  formula  developed  in  Ref.  1  for  the  headglow  radiant 
power, 

P  -  [0]  mN0  g(RU  (iN0^)  (Poo)-'4 

becomes  in  the  present  case,  where  T0j  =  const,  mj\fo  ~  Poj  and  g(M00) 
is  a  Mach  number  function  which  varies  little  for  the  present  set  of  data, 

p  -  40j)4  (p S'* 

Such  a  correlation  is  attempted  in  Fig.  lie.  The  deviation  from  the 
mean  is  larger  than  in  the  previous  figures.  This  is  partly  because  the 
formula  developed  in  Ref.  1  is  based  on  an  aerodynamic  model  substan¬ 
tially  different  from  that  obtained  in  the  test. 


7.2  CORRELATION  LIGHT  OUTPUT -MODEL  B 

A  typical  photograph  of  the  headglow  of  Model  B  is  shown  in  Fig.  12. 
The  sketch  in  Fig.  13  indicates  the  general  flow  characteristics.  A 
spherical  shell-like  area  around  the  contact  surface  between  the  two  gas 
streams  is  the  reaction  zone  and  produces  light. 
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The  results  of  part  of  the  runs  with  Model  B  are  tabulated  in 
Table  II  and  plotted  in  Fig.  14.  Figure  14  shows  the  effect  of  increas¬ 
ing  the  concentration  of  NO  in  the  model  jet  when  the  headglow  dimen¬ 
sions  are  kept  constant,  i.  e.  ,  when  P0j  and  POc0  are  kept  constant. 

Initially  the  total  light  output  is  proportional  to  m-^Q  as  would  be  ex¬ 
pected  for  Eq.  (2).  Note  that  riiQg  and  therefore  [0]  was  kept  constant 
for  the  runs  given  in  Fig.  14.  Further  increase  of  m-j^Q  does  not  in¬ 
crease  the  light  output;  a  "saturation"  point  has  been  reached. 

Only  two  runs  at  higher  T0  are  available,  apparently  in  the 
"unsaturated"  region  of  NO  mass  flow  (Fig.  14). 


SECTION  VIII 

EXPERIMENTAL  RATE  CONSTANT  FOR  HEADGLOW 


Several  photographs  taken  of  the  headglow  of  Model  B  were  analyzed 
to  determine  the  photon  production  per  unit  volume  in  the  headglow.  The 
exposure  of  the  film  is  related  to  the  number  of  photons  escaping  per 
unit  time  and  per  unit  area  of  the  projected  headglow.  By  scanning  the 
film  negative  as  outlined  in  Appendix  V  and  having  the  film  calibrated 
with  a  known  light  source,  the  photon  flux  per  unit  area  can  be  deter¬ 
mined.  Since  the  headglow  is  axially  symmetric  this  surface  flux  is 
related  to  the  radial  intensity.  The  inversion  from  surface  brightness 
to  radial  intensity  is  described  in  Appendix  VI.  In  calculating  the  radia¬ 
tion  intensity,  <E>  ,  the  measured  value  was  multiplied  by  2.  5,  accounting 
for  the  fact  that  the  film  response  curve  dropped  off  sharply  at  a  wave¬ 
length  of  0.  68  M.  The  factor  2.  5  was  obtained  by  comparing  the  film 
response  curve  with  the  spectrum  published  by  Fontijn  et  al.  (Ref.  3 
and  Fig.  1).  The  stagnation  point  values  of  4>  were  used  to  calculate  a 
rate  constant  according  to  Eq.  (1).  Half  the  value  of  the  concentrations 
[01  and  [NO]  on  either  side  of  the  contact  surface  were  used  in  calculating 
Ki  as  a  first  approximation  to  the  actual  concentrations  at  the  point  of 
maximum  light  intensity.  Therefore, 


l/2[NO]  l/2[o]  particle -sec 

For  the  range  [NO] /[0]  =  3.  6  to  54,  corresponding  to  the  first  linear 
part  in  Fig.  14,  and  with  [0]  =  1.  47  x  1012  atoms/cm3,  this  rate  con¬ 
stant  was  Kq  =  2.5  x  10“  12  cm3 /particle-sec  for  T0a>  =  290°K  and 
Po  =  50  m  Hg.  For  TOc0  =  800°K  and  Pq  =  40  M  Hg,  Kq  =  1  x  10-11  cm3 / 
particle-sec. 
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The  room  temperature  values  of  this  reaction  obtained  in  chemical 
flow  tubes  are  five  orders  of  magnitude  lower,  namely,  2,5  x  10“ 

(Ref.  7)  and  6.  4  x  10“-^  (Ref.  3). 

The  values  of  Ki  estimated  in  Ref.  1  from  free-flight  NO  release 
tests  were  from  1.1  x  10" 12  to  2. 6  x  10"  12  cm  3/  particle  -sec.  In  that 
case,  the  value  of  fo]  was  assumed  to  be  that  obtained  by  mass  spec¬ 
trometer  flights,  and  with  the  measured  light  output  an  effective  value 
of  Ki  was  calculated.  It  is  concluded  that  the  present  wind  tunnel  tests 
resolve,  at  least  within  an  order  of  magnitude,  the  problem  posed  by 
the  NO  release  flights  reported  in  Ref.  1.* 


SECTION  IX 

CONCLUDING  REMARKS 


The  present  test  series  indicated  a  much  higher  light  output  of  the 
reaction  of  NO  and  O  in  the  headglow  of  a  sounding  rocket  than  would  be 
expected  from  previous  laboratory  tests.  The  cause  of  this  discrepancy 
has  not  yet  been  determined.  Points  under  study  are  the  applicability 
of  laboratory  rate  constants  to  an  aerodynamic  flow  pattern  as  is  of, 
interest  here  and  the  differences  between  spectra  obtained  with  the  head- 
glow  and  laboratory  flow  tubes. 

A  method  of  producing  O  has  been  adapted  to  a  low  density  wind 
tunnel  flow  at  temperatures  well  below  the  dissociation  temperature. 

A  method  of  measuring  O  concentrations  in  a  low  density  super¬ 
sonic  flow  has  been  developed. 
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Fig.  1  Spectrum  of  Nitric  Oxide  Glow 
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Microwave  Cavities  Mounted  on  Oxygen  Intake  Manifold 
Fig.  3  Microwave  Antenna  and  Cavities 
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Fig.  4  Schematic  of  Low  Density  Tunnel  (Aerospace  Research  Chamber  (8V)) 
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Fig.  5  Sampling  Tube  Calibration  Arrangement 


AEDC-TR-66-105 


Chamber 


NO  Supply  Line  Wrapped  with 
Heater  Tape  and  Aluminum  Foil 

, _  tiplier  Tube 

[RCA  6217)  Mounted 

Note:  Entire  Tube  Wrapped  with 

Heater  Tape  and  Aluminum  Foil 


Except  at  Photomultiplier  Location 


Steel  Flange  Open 
to  Chamber 


Teflon  Tube  with  Small  Holes 


L  Teflon®!  ntake  Section 


Dimensions  in  cm 


Fig.  6  Sampling  Probe  Installation  -  ARC  (8V) 
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Fig.  9  Geometry  of  Test  Models 
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Fig.  11  Headglow  Light  Output,  Model  A 
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Symbols  as  in  Fig. 


Photomultiplier  Reading,  arbitrary  units 
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e.  Number  -  a(3(P0.) 


6  8  10 
5/4  gdvy 


(Poo^ 


20 

1/4 


'oo 


or~  [°Ja  NNO^coKmVi) 


1/4 


40  60  80 


3/4 


CO' 


(Eq.  13,  Ref.  1) 


Fig.  11  Concluded 
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Jet  Boundary  Shock 


Pressure  Distribution  along  Centerline 
(Not  to  Scale) 


Fig.  13  Sketch  of  the  Headglow  FI  ow  Field,  Model  B 
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Fig.  14  Headglow  Light  Output,  Model  B,  as  a  Function  of  NO  Flow  Rate  for 
Constant  Headglow  Geometry  (p  =  constant) 
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TEST 


Run  a. 
No.  % 


0, 

% 


Re /cm 


cm 


39 

40 

41 

42 

45 
44 

43 
57 

55 

56 

46 

47 

48 
54 

50 

49 
53 


0.50 


0.  50 


10.  55 

3. 

25 

3.  3 

1.  46 

10.  55 

3.  3 

1.  46 

10.  55 

3.  3 

1.  46 

10.  45 

3.4 

1.  42 

10.  45 

3.4 

1.42 

10.  45 

3.4 

1.  42 

7.  90 

3. 

31 

4.  3 

1.  13 

7.  90 

4.  3 

1.  12 

7.  90 

4.  3 

1.  12 

7.  95 

4.4 

1.  10 

5.  66 

3. 

39 

5.  9 

0.  85 

5.  66 

3. 

39 

5.  9 

0.  85 

5.  69 

3. 

39 

5.  9 

0.  85 

All  runs:  Tunnel  flow  T0  =  =  800°K 

WCD  Uq) 

Nitric  oxide  T0-  =  Tq.  =  290°K 
J  J 


I 

MODEL  A 


T 

X  CD’ 

°K 

P 

H  Hg 

PO’ 

M  Hg 

Poi 

p-^-x  1(T4 

1,  Arbitrary 
Units 

272 

1.42 

18.  5 

1.  81 

6.45 

1.  40 

18.  2 

6.  23 

36 

1.  40 

18.  2 

12.  0 

65 

1.40 

18.  2 

16.  8 

68 

25  7 

2.  82 

40.  0 

1.  71 

8.5 

257 

2.  82 

40.  0 

3.  50 

16.  75 

257 

2.  82 

40.  0 

6.  90 

27 

255 

2.  90 

41.  0 

6.  92 

35.  5 

255 

2.  90 

41.  0 

9.  75 

45 

255 

2.  90 

41.  0 

12.  0 

55 

250 

3.  55 

51.  8 

1.  24 

4.  1 

3.  60 

52.  1 

2.  46 

11.  55 

3.  60 

52.  1 

4.  78 

18.  75 

3.  60 

52.  8 

7.  20 

35 

242 

4.  50 

69.  0 

1.  93 

8,  25 

242 

4.50 

69.  0 

3.  54 

16.  25 

242 

4.  50 

69.  0 

5.  26 
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TABLE  II 

TEST  RESULTS,  MODEL  B 


Run  M  R_/nrn  T°°°’  T“’  P°>J  pV  mNO’  l>  Arbitrary 

No.  %  %  00  '  cm  °K  °K  m  Hg  M  Hg  mg/sec  Units 


7 

1. 

50 

5. 

84 

3. 

40 

19.  5 

0.  26 

2! 

30 

8 

8 

3. 

4 

52 

.  9 

77.  4 

39 

8 

6. 

05 

3. 

39 

18 

.  5 

0.  27 

3. 

3 

50.  8 

62.  0 

58 

9 

5. 

75 

19 

.  0 

0.  26 

3. 

35 

51 

.  8 

42.  2 

55.  5 

10 

' 

5. 

75 

18 

.  5 

0. 

27 

3. 

3 

50 

.  8 

25.  8 

54 

32 

] 

1. 

42 

6. 

13 

11.  0 

48 

25 

1. 

38 

6. 

13 

11.  0 

43 

26 

1.  42 

6. 

13 

8.  59 

53 

13 

1.  50 

5. 

75 

19.  0 

0. 

26 

3. 

35 

51 

.  8 

8.  60 

51 

23 

1. 

38 

6. 

13 

3. 

38 

17.  8 

0. 

28 

3. 

15 
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1. 

50 

5. 

75 

3. 

39 
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0. 

26 

3. 

35 

51.  8 

6.  38 

44 

27 

1. 

42 

6. 

13 

3. 

39 

18.  5 

0.  27 

3. 

3 

50.  8 

6.  12 

42 

24 

1. 

38 

3. 

38 

17 

.  8 

0. 

28 

3. 

15 

48 

.  5 

6.  03 

56.  5 

28 

1. 

42 

3. 

39 

18.  5 

0. 

27 

3. 

3 

50.  8 

3.  68 

25 

29 

1. 

42 

3. 

39 

18 

.  5 

0.  27 

3. 

3 

50 

.  8 

1.  13 

5.  75 

30 

1. 

42 

3. 

39 

18 

.  5 

0.  27 

3. 

3 

50.  8 

0.  245 

3.  75 

33 

0.  50 

9. 

85 

3. 

25 

3 

.  3 

1. 

46 

800 

25  7 

2. 

85 

40.  0 

11.  0 

13 

34 

0.  50 

9. 

85 

3. 

26 

3 

.5 

1. 

38 
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2. 
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42 
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3.  66 

4.  8 

All 

runs 

:  T0- 
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T  ' 
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APPENDIX  I 

DISSOCIATION  OF  NITROGEN  TETROXIDE 


For  the  general  case 


2 

i  =  1 


vi  Ai: 


m 

2 

i  =  1 


vi"  Aj 


the  law  of  mass  action  gives  for  the  equilibrium  constant 

K«  .  4A  -  K«(X)  -  -fill - - 


K 


R 


-g-  (Pi')  ’'i 


Thus  for  the  reaction 


we  have 


Kd 

2  N02  2=TT  N204 


Ke(T) 


(  Pno2  ) 2 

(pN204  ) 


Also  when  expressed  in  the  degree  of  dissociation 

=  [no2] 

[no2]  +  2  [n2o4] 

there  results 

4a*  Kb(T) 

"  P 


where  P  is  the  total  pressure. 
Equation  (1-2)  can  also  be  written  as 


an  =  an(P’  T) 


‘4 


1  +  4 


Ke(T)  J 


(1-1) 


(1-2) 


(1-3) 


Note  that  an  is  a  function  of  P  and  T  but  that  Ke  is  a  function  of  T  only. 

Equilibrium  composition  data  for  reaction  Eq.  (1-1)  are  collected 
in  Ref.  1-1.  Since  low  pressure  (P  «  1  atm)  data  were  not  readily 
available,  the  values  of  Ke  were  calculated  from  the  tables  in  Ref.  1-1, 
using  Eq.  (1-2)  and  plotted  in  Fig.  1-1.  The  straight  line  correlation 
in  Fig.  1-1  gives 

Ke  =  3.01  x  109  exp  (1-4) 


39 


AEDC-TR-66-105 


This  value  was  used  in  Fig.  1-2  to  calculate  an  =  a;n(P,  T)  using  Eq.  (1-3). 
Finally  the  degree  of  dissociation  is  given  as  a  function  of  pressure  at 
room  temperature  in  Fig.  1-3. 


REFERENCE 


1-1  Nitrogen  Tetroxide.  Allied  Chemicals  Co. 


40 


o  P  =  0. 5  atm  Obtained  from 
O  1  Equilibrium  Data 


K  in  Dof  1.1 
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APPENDIX  II 

CALCULATION  OF  MASS  FLOW  OUT  OF  FLASK 


The  mass  flow  out  of  the  flask  was  calculated  from  a  pressure  time 
recording.  The  relation  between  m  and  —  is  derived  here. 


The  mass  flow  is 


(II- 1) 


where  V  is  the  volume  of  the  flask. 

For  dissociated  gas, 

p-PFT(1  +  %>  (II- 2) 

where  M  is  the  molecular  weight  of  the  undissociated  gas. 


Extrapolating  the  relaxation  times  listed  in  Ref.  1-1  and  plotted  in 
Fig.  II- 1,  the  relaxation  time  for  dissociation  is  estimated  between  10"^ 
and  10~6  sec  for  the  pressure  range  of  interest  for  NO2  (P  =  4  to  100  mm  Hg). 
Therefore,  the  equilibrium  degree  of  dissociation  of  Appendix  I  is  assumed. 


Logarithmic  differentiation  of  (II -2)  gives 


(r)©-(r)e-)-(T^)(^j 


The  term 


dan 

dt 


is  expressed  in  terms  of  P  and  T  with  Eq.  (1-3) 


(II- 3) 


dt 


and  since  according  to  Eq. 


(1-4) 

K  =  ae~b/T 

e 


and  finally 


44 


A  E  DC-T  R-66- 105 


Substitution  of  the  above  formulas  in  Eq.  (II- 1)  gives 


-  V 

dP 

VM  (/ 

a  a2  \ 

,  _ n_  _ _ 2_) 

dp 

dt~ 

RT  (1  +  an)  j  \ 

2  2  / 

dt 

p 

T 

1 

\ 

+  l/2an  (1  -  an)  — ] 

d’I’| 

TTj 

(II-4) 

For  the  system  N2O4  tz;  2NC>2  the  value  of  b  =  7120°K. 

For  a  typical  run  in  the  present  test,  At  -  20  sec,  P  =  40  mm  Hg, 
a  -  0.  5  at  T  =  300°K,  and 


dP  ~  ft  l  mm  Hg 
dt  sec 


The  two  terms  in  the  square  bracket  of  Eq.  (II-4)  become 

0.1125  -  0.53 

dt 


The  requirement  for  the  last  term  to  be  negligible  is  then 

il  «  0.21 

dt 


or  for  At  =  20  sec 


|  AT|  «  4.2°K 


This  condition  may  not  always  be  fulfilled  so  that  the  completed  Eq.  (II-4) 
may  have  to  be  used. 


In  laboratory  tests  a  thermometer  in  the  flask  showed  no  change  of 
temperature,  i.  e.  ,  AT  -  0.  Thus  for  NO2  in  the  present  test  the  mass 
flow  out  of  the  flask  is  given  by: 


m(N  02)j 


VM  2  -a  dp 

_ _ n _  _ 

RT  2  dt 


(II -5) 


Using  V 
T 
M 


5000  cm^ 

293°K 

92  gm/mole  (N2O4) 

rn(N02)j  =  ~  25,2 


dP  mg 

dt  9  sec 


where  P  is  in  mm  Hg 


In  the  case  of  NO  (a  =  0),  M  =  30  gm/mole  the  mass  flow  out  of  the 
same  flask  is 


mN  0 


-8.21 


mg 

sec 
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APPENDIX  III 

MASS  FLOW  THROUGH  SAMPLING  TUBE  AND  INTAKE  EFFICIENCY 


The  pressure  drop  along  the  sampling  tube  is: 


dp 

dx 


A 

d 


1 

2 


pn 


2 


and  assuming  laminar,  Poiseuille  flow  the  friction  factor  is 

^  _  64  _  64p 

Re  p  ud 


In  terms  of  the  mass  flow  mst 


so  that 


dp 

dx 


P2  -  P2  = 
1  2  1  1 


pu  JL.  d2  through  the  tube  this  becomes 
128  mRT 

77  ^  PMd4 


256  . 

-  mst 

17 


RTp 


(xi 


X2) 


or 


(III-l) 


The  validity  of  Eq.  (III-l)  was  checked  with  the  sampling  tube  mounted 
outside  the  chamber  and  gas  supplied  by  a  bottle.  The  mass  flow  was 
then  measured  with  a  mass  flowmeter  (Fig.  5).  The  value  of  mst 
calculated  from  Eq.  (III-l)  agreed  within  10  percent  of  that  obtained 
with  the  flowmeter. 


o 

In  the  tunnel  runs,  the  ratio  8  =  - - - is  used.  The  molecular 

mQ  +  mN 

weight  M  is  then  (/3  4  1) 

mO  M0 

M  =  — —^—7—^ - ;  2  2-  =  /3Mn  +  (1  -  p)  Mn  =  4(7  +  fi) 

mO,  +  mN,  2  2 

Furthermore,  during  the  tunnel  runs  Tw  =  290°K  for  both  heated  and 
unheated  runs.  And  using  present  tube  dimensions,  the  final  mass  flow 
formula  is 

ihst  =  0.68  x  10-3  (7  +  /3 )  (Pi  -  P2)  mg/sec  (III -  2 ) 

where  Pj_  and  P2  are  in  m  Hg. 
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The  mass  flow  through  the  tube  as  calculated  from  the  measured 
pressures  with  Eq.  (Ill - 2)  is  compared  with  the  maximum  mass  flow 
( u^ A intake )  in  Fig.  Ill  - 1 .  The  plot  of  efficiency,  p,  versus  Knudsen 
number  (mean  free  path  divided  by  tube  intake  diameter)  shows  that  the 
tube  is  operating  in  the  transition  regime.  At  low  Kn,  continuum  flow, 
the  shock  there  is  apparently  an  attached  normal  shock  and  p  approaches 
100  percent.  The  limit  for  Kn  »  1  is  not  known,  but  the  trend  shown 
in  Fig.  Ill  - 1  is  reasonable.  For  Kn  >  0.  1  the  shock  becomes  very 
thick,  and  substantial  spilling  will  occur.  The  extent  of  the  spilling  will 
depend  on  the  conditions  downstream  in  the  tube. 

For  the  higher  temperature  runs  (T0  =  800°K  and  (Tw/T0)  -  0.  36) 
the  intake  efficiency  is  increased  at  a  given  Knudsen  number  because 
of  the  effect  of  wall  cooling  of  the  gas  flowing  through  the  tube  or  at  the 
inlet. 
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APPENDIX  IV 

WALL  RECOMBINATION  IN  A  GLASS  TUBE 


is: 


The  general  scheme  for  recombination  of  O  proposed  in  Ref.  IV- 1 


0  +  wall  >  1/2  02  +  wall 

0  +  02  +  M  03  +  M 


(IV-1) 

(IV-2) 


0  +  03  2  02  (IV- 3) 

The  number  of  O  atoms  striking  the  wall  per  unit  time  in  a  section  dx  is 


(1/4  N  ua)  area  =  ~  N  ua  rdx 


where  N  is  the  number  of  atoms  per  unit  volume  and  ua  is  the  arithmetic 
average  atomic  velocity.  Of  these,  the  fraction  7  recombines  at  the  wall 
according  to  Eq.  (IV-1).  Expressing  this  in  standard  reaction  rate  nota¬ 
tion,  we  have 

— —  N  u  rdx  =  K,  N  ttt2  dx 
2  a 


or 


(IV-4) 


Note  that  Kq  has  the  dimension  [T"l]. 

The  following  values  of  7  have  been  reported  in  the  literature. 


Elias,  Ogryzlo,  and  Schiff 

7  =  0.  77 

x  IO-4 

Linnett  and  Marsden 

7=1.2 

x  10"4 

Kaufman 

7  =  0.  2 

x  10"4 

Kistiakowsky  and  Volpi 

7  =  0.  2  x  10  4  (deduced 
minimum  value) 

Herron  and  Schiff 

7=1.1  x  10“4  (deduced 
upper  limit) 

Obviously  the  value  of  7  will  depend  on  the  condition  of  the  wall.  So- 
called  poisoning  of  the  wall  with  a  suitable  chemical  will  reduce  wall 
recombination. 
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From  Eqs.  (IV-1),  (IV-2),  and  (IV-3) 


d[o] 

at 


=  K,  [0]  +  K2  [0]  [OJ  [M]  +  K,[0]  [OJ 


JoJ 

d  t 


=  ~Ka  [0]  [OJ  [M]  +  K3  [0] 


[OJ 


(IV-5) 


It  is  now  assumed,  as  in  Ref.  IV-1,  that  the  ozone  concentration  is  con¬ 
stant,  thus 

K3  [0]  [OJ  =  K2  [0]  [OJ  [M] 


Substitution  in  Eq.  (IV-5)  gives 

d  in  [o] 


dt 


or 


-  In 


[  [0]  1 

1 

X 

II 

,o 

=  K,  +  2K2  [OJ  [M] 


K,  +  2Kj  [OJ  [M] 


(IV-6) 


where  a  uniform  velocity  u  is  assumed  along  the  tube 


Since  the  gas  in  the  tube  was  at  room  temperature. 


[OJ  - 


(P)o, 

760 


L0,  particle/cm3 


(P) 

[M]  =  ~~~  L0,  particle/cm3 


Now  for  the  present  tests, 

<p>o. 


<P>r 


111  n 

_  _ — =  /3 

'M  mO,  +  mM 


so  that 


[OJ  DM]  =  /3  (1  -  /3)  Li 


(7  60) 


(IV- 7) 


In  the  three-body  Eq.  (IV- 2)  the  third  body  M  can  be  also  O2  so  that  for 
j3  =  1  the  product  [OJ  [W]  should  be  [OJ2 .  In  the  present  test  for  most 
cases:  j3  «  1  so  that  /3(1  -  /3)  «  j3  and 

[OJ  EM]  -  13  (P 

The  following  numerical  values  were  used: 

y  =  10"4 

ua  =  4.44  x  104  cm/sec 
r  =  3.5  cm 

L0  =  2.687  x  1019  particle/cm3 


K, 


0.634  sec  ’ 
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From  Ref.  IV- 1, 

K2  =  10 14  -J-  =  2.8  x  10-34  - ^ - 

mole  sec  particle  -sec 

so  that  Eq.  (IV-6)  becomes 

- In  -4^3 -  =  (0.634  +  0.70/3P2) 

t°lx  =  o 


where  P  =  the  pressure  in  mm  Hg 


(IV-8) 


Equation  (IV-8)  is  plotted  in  Fig.  IV- 1.  This  figure  shows  that  for 
average  flow  velocities  >10  m/sec  the  recombination  ifcfffe  tube  is  less 
than  a  few  percent. 


Note  that  the  difference  between  the  curves  P — >0  and  P  4  0  repre¬ 
sents  the  effect  of  volume  recombination. 
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Fig.  IV-1  Oxygen  Recombination  in  a  Glass  Tube  with  Uniform  Velocity 
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APPENDIX  V 
PHOTODENSITOMETER 


The  photographic  density,  D,  of  a  negative  is  defined  in  terms  of 
the  intensity  of  a  light  beam  transmitted  through  a  negative: 

D  =  log  -is- 

where  I  =  intensity  of  transmitted  light  beam 
IQ  =  intensity  of  undisturbed  light  beam 

The  film  characteristic  is  a  curve  of  D  versus  log  E,  where  E  is  the 
exposure  (in  ergs  or  photons)  of  the  negative.  This  curve  is  character  ¬ 
istic  for  the  film  used.  In  the  present 
investigation  each  film  roll  was  calibrated 
by  exposing  the  film  by  a  series  of  known 
exposures,  using  a  light  box  which  was 
calibrated  against  a  standard  source  at 
AFCRL.  The  light  source  was  calibrated 
at  various  frequencies.  Its  spectral  distri¬ 
bution  matched  the  film  frequency  character¬ 
istic.  During  this  film  calibration,  the 
light  box  was  placed  at  the  same  distance 
from  the  camera  as  the  distance  from  model  to  camera,  so  that  no  geo¬ 
metric  adjustment  was  required. 


The  negatives  were  run  through  a  densitometer  as  sketched  below. 
The  intensity  of  the  transmitted  light  beam,  along  the  dotted  line,  was 
recorded  on  a  strip  chart  recorder. 


□  □□!□□□□□ 


The  following  readings  are  noted  (Ref.  V-l): 

R0  =  full  opaque  reading  of  densitometer 
Rb  =  reading  outside  headglow,  background  reading 
Rm  =  reading  in  the  headglow,  measurement 
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I  R,  -  R 

D  =  log  -j—  =  log  ^ 


-  R 


(V-l) 


In  the  linear  range  of  the  film  characteristic, 

=  Pi  ~  P2 
y  log  E,  -  log  E2 


(V-2) 


and  substituting  (Ej)  in  (E2), 


(V-3) 


Either  1  or  2  can  be  taken  as  a  reference  condition 

Em  _  /  Rref  ~  Hq  \ 

Eref  \  Rm  -  R0  / 

In  the  present  case  Rref  =  Rb  was  taken,  i.  e. ,  the  background  level  of 
the  film. 


Era  1 

e„  '  TW 


Efo  is  then  immediately  read  of  the 
film  characteristic  and  the  final 
formula  is 


1  < 


Em 

Eb 


<  OO 


(V-4) 


r 


Rm  -  R0  /  / 

— - -  ,  0  <  r  <  1 

Rb  -  R0  ~  “ 
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APPENDIX  VI 

RADIAL  INTENSITY  DISTRIBUTION  DERIVED  FROM  PHOTOGRAPHS 


For  an  axially  symmetric  geometry  the  surface  brightness 
(erg/cm^-sec)  as  shown  in  the  film  exposure,  can  be  related  to  the 
specific  radial  intensity  distribution  (erg/cm^-sec).  For  the  present 
case  the  following  assumptions  are  made: 

(a)  The  gas  in  the  headglow  is  nonabsorbing  and  nonscattering 

(b)  The  gas  is  either  a  grey  medium  or  all  parts  of  the  headglow 
emit  in  the  same  frequency  range 

(c)  The  distance,  d,  from  the  film  to  the  headglow  is  such  that 

d  »  R,  so  that  a  point  on  the  film  receives  radiation  from  one 
column  of  headglow  parallel  to  the  y-axis  only. 


For  the  headglow  the  following  coordinate  system  is  introduced: 


x 


y 


I  (z, 


y 


X 


If  (x,  C  )  is  the  projected 
surface  brightness  of  a 
dx  column  at  x  =  constant, 

and  I  (z,  r)  is  the  radiation 
emitted  per  cm^  per  second 


then 

C(x,  z)  =  /  I(z,  r)  dy 


or,  changing  variables  with  x2  +  y2  =  r2, 

{(»■■)•(  r‘<a'  'V'  (VI-1) 

Vr  *-  x 

where 
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The  limits  of  integration  are  r  =  R  and  r  =  x  and  this  range  is  covered 

Em 

twice.  From  the  film  the  quantity  F(x,  z)  =  is  determined 

(Appendix  V)  which  is  proportional  to  C  (x,  z).  If  all  distances  are  non- 
dimensionalized,  Eq.  (VI- 1)  becomes  at  a  given  value  of  z, 


r.  ,  s  „  ,R  rl  (r)  dr 

F  (x)  =  2  f  - - 

X  V  r  -  X 


(VI-2) 


In  the  present  case,  the  reference  length  was  the  model  diameter,  DmodeL 
read  from  the  photograph,  so  that  in  the  notation  of  Eq.  (VI- 2)  the  radia¬ 
tion  emitted  per  cm3  per  second  is 

I(r)  Ej, 


o  d  el 


Now  for  a  given  distribution  F(x)  at  a  fixed  value  of  z  the  distribution 
I(r)  is  desired.  It  is  shown  in  Ref.  V-l  that  the  Abel  integral  equation 


f  (s)  =  / 

a  (  s  -  t  )a 


dt 


has  the  inversion 


{  (s) 


sin  an  d  -l 

7T  d  t  l  (t-s)1_a 


ds 


T  (  \  d  fR  F  (x)  dx 
77  r  1  (r)  =  -  J 


After  substitution  of  x2  — >s  and  r2— *t  with  a  =  1/2,  Eq.  (VI-1)  becomes 

, _  (VI-3) 

Vx2  -  r2 

This  is  also  shown  in  Ref.  VI-2.  The  last  equation  is  now  written  in  a 
different  form: 

R  /  \  /  p  x=R 

TrrI  (r)  =  /  F  (x)  d  (Vx2  -  r2  J  =  |  [_F  (x)  Vx2  -  r2J 

-  /  Vx2  -  r2  dF(x)|  =  |f(R)  V'R2  -  r2  -  /  Vx2  -  r2  F'(x)dxj 

In  the  present  case  F(R)  -  0  so  that 


j  z  \  d  R  .  /  j  2  t-  -  /  \  j  R  r  F '  (x)  dx 

77  r  I  (r)  =  -  J  Vx  -  r  F  (x)dx  =  f  - — 


dr  r 


V  X2  - 


and  finally. 


I  1 

I(r)  =  -±-  f 


F  '  (x)  dx 


„  r  V  X2  -  r2 


(VI-4) 
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Since  the  integral  becomes  infinite  at  the  lower  integration  limit 
this  form  is  not  convenient  for  numerical  calculations.  The  problem  is 
circumvented  by  writing 


Now 


I(r)  = 


1 

-  J 

77  r 


R 


F'(x)  -  FTP  dx  +  J_  rR  F  '  (r)  dr 
\/x2  -  r2  77  r  Vx2  -  r2 


R 

/ 


r 


x  +  V  x2  -  r2 


for  x2  > 


The  latter  condition  is  satisfied  since  in  our  case  R  >  x  >  r. 


The  final  result  is: 


I  (r)  = 


F  '  (x)  -  F  hr) 


dx 


F  '  (r) 


In 


Vx2 


R  +  Vr2 


(VI-5) 


In  general  the  accuracy  of  numerical  calculation  is  increased  if 
derivatives  of  curves  through  experimental  data,  F(x)  in  this  case,  can 
be  avoided.  Returning  to  Eq.  (VI-4), 


77 1(r)  =  f 

r 


■H  d[F  (x)] 


[F(x)  -  F (r)] 


F  (R)  -  F  (r) 

Vr2  ~  r2" 


F  (x)  -  F  (r) 


v: 


+  / 
r 


[F  (x)  -  F  (r)  ]  xdx 
(x2  -r2)3/2 


(VI-6) 


Assuming  for  the  moment  that  the  second  term  is  zero,  the  result  is 


I  (r)  = 


1  F(R)  -  F  (r) 

TT 


[f  (x)  —  F  (r)]  xd  x 
(x2  -  r2)3/2 


(VI-7) 


The  second  term  in  Eq.  (VI-6)  can  be  written  as 


F  (x)  -  F  (r) 

=  lim 

F  (r  +  f)  -  F  (r) 

=  lim 

F  '  (r  +  e)  -  F  '  (r) 

1/2  [2  (r  +  e)  -  2r] 

V?  -r2 

x  =  r  e->  0 

V.(r  +  e)2  -  r2 

£-*  0 

V(r  +  6)2^7 


•=  lim  JL'b  +  t) 

£-»  0  f 


V  (r  +  £)2  -  r2 


But  since  t^dr  the  first  factor  is  also 

L  =  lim  F "  (r)  Vt2  +  2  V  =  0 
£-*  0 


AF  (r)  p”  (r) 
Ar 


so  that 
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